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ABSTRACT 

The Large Observatory For X-ray Timing (LOFT), selected by ESA as one of the four Cosmic Vision M3 
candidate missions to undergo an assessment phase, will revolutionize the study of compact objects in our 
galaxy and of the brightest supermassive black holes in active galactic nuclei. The Large Area Detector (LAD), 
carrying an unprecedented effective area of 10 m 2 , is complemented by a coded-mask Wide Field Monitor, in 
charge of monitoring a large fraction of the sky potentially accessible to the LAD, to provide the history and 
context for the sources observed by LAD and to trigger its observations on their most interesting and extreme 
states. In this paper we present detailed simulations of the imaging capabilities of the Silicon Drift Detectors 
developed for the LOFT Wide Field Monitor detection plane. The simulations explore a large parameter space for 
both the detector design and the environmental conditions, allowing us to optimize the detector characteristics 
and demonstrating the X-ray imaging performance of the large-area SDDs in the 2-50 keV energy band. 
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1. INTRODUCTION 

The LOFTjll WFM is a coded aperture imaging experiment designed on the heritage of the SuperAGILE ex- 
periment P successfully operating in orbit since 2007.3 With the ~120 /itm position resolution provided by its 
single-sided Silicon microstrip detector, SuperAGILE demonstrated the feasibility of a compact, large-area, light, 
low-power and high resolution X-ray imager, with steradian-wide field of view. The LOFT WFM applies the 
same concept, with improvements provided by the superior performance Silicon Drift Detectors (SDDs) in place 
of the Si microstrips. These detectors provide a lower energy threshold and better energy resolution with respect 
to the microstrips, and represent the state-of-the art in large area, monolithic Silicon detectors. 

The working principle of the WFM is the classical sky encoding by coded masks^ and is widely used in space 
borne instruments (e.g. INTEGRAL, RXTE/ASM, Swift/BAT). The mask shadow recorded by the position- 
sensitive detector can be deconvolved by using the procedure^ to recover the image of the sky, with an angular 
resolution that depends on the mask element pitch, the mask-detector distance and the detector spatial resolution. 
In this paper we describe the results obtained by means of a Monte Carlo simulator specifically developed to 
evaluate the spatial resolution capabilities of the large area Silicon Drift Detectors. This tool allowed us to 
understand and optimize the detector performance by investigating different detector designs (i.e. anode pitch), 
environmental conditions (detector temperature, radiation damage) and read-out electronics performance. 
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In the following we will describe the SDD working principle and the charge cloud dynamics inside the Silicon 
bulk (Section [2]), the preliminary imaging results obtained with a 8-channels discrete read-out (Section [3]), the 
structure of the MC simulator developed for the LOFT/WFM (Section]!]) and the simulation results (Section|5|. 
Finally, in Section [6] we draw our conclusions. 




2. WORKING PRINCIPLE OF THE SILICON DRIFT DETECTOR 

As reported above, the proposed technology for the WFM detection plane is represented by an improved version of 
the large-area Silicon Drift Detectors (SDDs) developed for the Inner Tracking System in the ALICE experiment 
of the Large Hadron Collider (LHC) at CERNJM3 by one of the scientific institutes in the LOFT Consortium — 
INFN Trieste, Italy — in co-operation with Canberra Inc. The key properties of the Si drift detectors^ are their 
capability to read-out a large photon collecting area with a small set of low-capacitance (thus low-noise) anodes 
and their very small weight (~1 kg m~ 2 ). 

The working principle of the Silicon drift detectors is shown in Figure [l] When a photon is absorbed in 
the Si bulk it creates a cloud of electron-hole pairs. The holes are collected on the cathodes implanted on both 
the detector surface while the cloud of electrons is focused in the middle plane of the detector by a parabolic 
potential distribution. The e~ cloud then drifts towards the read-out anodes by means of a constant electric 
field sustained by negative voltages applied to the cathodes with progressively reducing amplitude down to the 
anodes at V. When focused in the detector middle plane, the cloud quickly assumes a Gaussian shapeP The 
diffusion in Si causes the electron cloud to expand by a factor depending on the temperature, the electric field 
and drift length, as described by the relation: 



a = ^2Dt + al = f^f^ + °t = + -I « 

where D is the diffusion coefficient, t is the drift time, kg is the Boltzmann's constant, fj, is the electron mobility, 
T is the absolute temperature, uo is the initial cloud dimension and x is the drift length. In Equation [T] we used 
the Einstein relation D = /iksT/q to derive the final result. 

The charge distribution over the collecting anodes then depends on the (x,y) absorption point in the detector. 
In fact the z-th anode collects a charge Qi which is the integral of the Gaussian cloud between the anode 
boundaries. Qi can thus be expressed as: 
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where K is a value, changing event by event, that represents a common baseline (i.e. a common mode noise, see 
|10| ), Qtot is the total cloud charge, ]ji is the center of the i-th anode, yo is the photon absorption position in the 
anodic direction, p is the anode pitch, a represents the standard deviation of the charge distribution and Erf is 
the error function. 

The design of the large-area SDD is illustrated in Figure fijb). The Si tile is electrically divided in two halves, 
with two series of read-out anodes at the edges and the highest voltage along its symmetry axis (indicated as 
a solid line in Figure [T]d). The drift length is 35 mm, and the drift field is 360 V cm' 1 , which corresponds to 
1300 V maximum voltage. The drift velocity is then about 5 mm /is -1 which translates in a maximum drift 
time of ~7 /is. The maximum size of the charge cloud reaching the anodes can be estimated using Equation [T] 
with x — 35 mm (corresponding to an event absorbed at the bottom of the drift channel). This leads to about 
1.3 mm (97% containment radius) at room temperature. 




Figure 2. ALICE detector integrated with the discrete front-end electronics. 

3. THE ALICE HERITAGE 

The large area multi-anode SDD was initially designed in the framework of the LHC- ALICE experiment at CERN, 
and it demonstrated a spatial resolution better than 30 /im in the two dimensions in localizing the impact point 
of ionizing particles!^ The ALICE-D4 SDD features a sensitive area of 53 cm 2 , with the two detector halves 
read-out by means of two arrays of 256 anodes with 294 /im pitch. These detectors were produced on 5-inch 
diameter, 300 /im thick NTD (Neutron Transmutation Doped) Si wafers. Since 2008, 260 SDDs (for a total area 
of 1.37 m 2 ) have successfully been in operation at ALICE. 

Over the last few years R&D work has been carried out to characterize and optimize the same detector design 
for detection of soft X-rays. The preliminary results obtained under the XDXL INFN project (in collaboration 
with INAF, Politecnico di Milano, Universit di Pavia) with a spare detector of ALICE with a bread-board read- 
out based on discrete electronics, showed high spectral performance and good position resolution already at room 
temperature. A preliminary evaluation of the spectroscopic and imaging performance of the ALICE Silicon Drift 
Detector performed at the INAF/IAPS X-ray facility^ can be found in [ToJ[l3j[l4] . 

Figure [2] shows the SDD+FEE setup used in the measurements at the INAF/IAPS laboratory. Eight con- 
tiguous anodes of the SDD were separately connected to the same number of low gate capacitance JFETs 
(Cgs = 0.4 pF) used as the input transistor of Amptek A250F-NF charge sensitive amplifiers. The feedback 
capacitor Cf — 50 fF and a reset transistor are both integrated on the JFET die, allowing to reduce the input 
stray capacitance for a better noise performance. Nonetheless, some stray capacitance is unavoidable due to the 
need to adapt the small SDD anode pitch to the much wider one of the electronics. Such a parasitic capacitance 



was estimated to be as high as ^2.5 pF with the ALICE-spare detector set-up. Figure IJa) shows the ALICE 
SDD (equipped with the discrete FEE) imaging capabilities at T room for a 4.5 keV monochromatic X-ray beam 
(with a spot size of ^100 /im FWHM) as measured during the laboratory tests. It is worth noticing that, besides 
the noise performance of the discrete read-out electronics, also the small number of anodes (8) collecting the 
charge cloud introduces some systematics in the determination of the Gaussian cloud parameters (Equation [2| , 
both with respect to the spectroscopic (see Figure 7 in 10 ) and the imaging performance. 




(a) (b) 

Figure 3. (a) Reconstructed detector image at T=293 K for a 1800 s integration with a 4.5 keV monochromatic X-ray 
beam of ~100 /im FWHM, at yo ~ and xo ~ 10 mm. (b) Simulation of a Gaussian beam with 100 ^m FWHM centered 
at yo = mm and xo = 10 mm. 



4. THE SDD SIMULATOR 

In order to estimate the spectroscopic and imaging capabilities of the Silicon Drift Detectors and to verify the 
consistence between the experimental results and the assumptions, we developed a Monte Carlo simulator written 
in IDL language and describing the charge drift and diffusion inside the SDD as well as the detector-FEE noise 



characteristics. A first version of the simulator, discussed in 10 13 14 , allowed to model the ALICE-D4 SDD 
equipped with the 8-channels discrete FEE. The simulation results showed a very good agreement with the 
preliminary experimental measurements performed at T Toom , demonstrating the accuracy of the SDD modeling. 
In Figure |3jb) we show, as an example, the simulator output for the 8-channels discrete read-out at room 
temperature. This simulation was carried out considering the same experimental conditions of the measurement 
showed in Figure |3][a) . 

A new version of the simulator was then developed to study the performance obtainable by using an ASIC 
based FEE and by optimizing the SDD design (especially in regard to the anode pitch). The estimation of 
the imaging (and spectroscopic) capabilities of the SDD was carried out by means of simulations performed at 
different system temperatures and anode pitches. For these simulations, we considered a 32-channel integrated 
read-out (as baselined for the LOFT/WFM ASIC) and a noise figure consistent with the post-layout simulations 
of a test ASIC being developed within the XDXL project. The following scheme provides a step-by-step description 
of the simulator: 



1. A photon of energy E is generated in a position Xq (along drift direction) and yo (along anode direction) 
in the SDD; 

2. The photon is absorbed in the Silicon bulk and generates an electron cloud in (xo,yo). The estimation of 
the total charge in the cloud takes into account the average Silicon electron-hole pair generation energy 
(3.6 eV) and the Fano factor (0.115); 
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Figure 4. (a) Number of anodes collecting the electron cloud (99.7% containment radius) as a function of the photon 
absorption point along the drift direction for anode pitches from 145 /im to 294 /im. The curves are simulated considering 
HV=-1300 V, T=253.15 K, cr =0 Mm. (b) Expected ENC at °C (red) and -20 °C (black) for a 145 /im pitch SDD 
in function of the shaping time. Curves for both Begin of Life (BOL, dashed line) and End of Life (EOL, solid line) 
performance are shown. 



3. The cloud, which can be assumed to have an initial width of <7o ~ 0, is propagated along the drift direction, 
reaching the read-out anodes with a width of c, which is a function of the diffusion coefficient D and of 
the drift time t as described in Equation [l] 

4. The electron cloud is then collected by a number n of anodes depending on the width of the cloud and 
on the photon conversion point (xQ,y ). The charge fraction Q t which reaches each anode is described 
by Equation [2] Figure Ha) shows the number of anodes collecting the electron cloud as a function of the 
photon absorption point along the drift direction for anode pitches from 145 jum to 294 /zm with T = 253 K 
and E = 360 V cm -1 . 

5. At this stage, the simulator takes into account the Equivalent Noise Charge (ENC) for each anode indepen- 
dently. The ENC is estimated considering an integrated FEE with the following characteristics: CR-RC 
shaper, Cf cc dback = 21 fF, C s t ra y = 200 fF, and an available power per channel of 0.722 mW. The anode 
capacitance is estimated from the anode pitch and the detector thickness (between ~ 50 fF and ~ 130 fF 
for pitches in the range 145-294 /im and 450 /zm thickness) while the leakage current is evaluated consid- 
ering the intrinsic current measured in the laboratory plus the contribution due to the radiation damage 
after 5 years in orbit at 600 km altitude and 5° inclination, estimated by means of simulations performed 
with the ESA's Space Environment Information System (SPENVIS, http://www.spenvis.oma.be). The 
dependence of the leakage current on the temperature is taken into account by considering the relation 

I icak (T)ocT 2 -e- E ^ 2k - T (3) 

where T is the absolute temperature, q the electron charge, E g the Si bandgap (1.2—1.3 eV in the temper- 
ature interval 250 < T < 295 K) and ks the Boltzmann's constant. In Figure Qb) we show the expected 
ENC at °C and -20 °C for a 145 zzm pitch as a function of the shaping time for both the Beginning Of 
Life and the End Of Life cases. A Common Mode Noise (CMN), with an rms value of 25 e~ is then added 
to the 32 ASIC channels. This CMN is due to the presence of noise on the detector bias line during the 
laboratory tests and can be reduced by proper filtering. 

6. The charge distribution is eventually fitted using Equation [2j The fitting parameters Qtot, yo an d cr (and 
K) allow to completely characterize the physical properties of the electron cloud (i.e. E, xq, yo). 



Table 1. Parameters used in the detector imaging simulations. 
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The simulations were performed exploring a large parameter space. Table [T] summarizes the different values 
used in the simulations concerning the anode pitch, the detector temperature, the incident photon energy and 
the photon absorption point in the anodic (y) and drift (x) coordinates. 

5. ANALYSIS OF SIMULATED DATA 

For each simulation configuration we generated a large number (2 x 10 5 ) of photons. The output is saved as 
FITS tables, including a truncation of the values to simulate the ADC response. The event list containing the 
simulation results were then analysed fitting Equation [2] by means of a least squares minimization algorithm 
written in C (CMPFIT, http://www.physics.wisc.edu/~craigm/idl/cmpfit.html). 

Figure [5](a) and Figure [5jb) show the anodic and drift spatial resolution (FWHM) for a 2 keV photon at 
253 K for different anode pitches, for several photon absorption points in the drift direction and for two different 
positions in the anodic direction: at the center of the anode (y=0), and 50 fim far from the anode center. The 
results shown in the plots indicate that a pitch smaller than 175 fim is required to minimize the systematic effects 
on the position reconstruction introduced by the discretization of the charge cloud read-out and thus to optimize 
the detector spatial resolution at low energies (<2. 5-3.0 keV). For energies larger than 3 keV, the Signal-to-Noise 
ratio in each bin of the integrated Gaussian cloud increases, thus improving the spatial resolution as shown in 
Figure [5](c) and Figure [5](d) for E ph =3.5 keV and in Figure [6] for E ph =6.0 keV. 

Looking at the simulation results in Figure [5] and Figure [6] is it clear that a small anode pitch translates in a 
better spatial resolution for photons absorbed in the first part (~ 10-15 mm) of the drift channel. On the other 
hand, the dependence of the anodic and drift resolution on the pitch is less important at the end of the drift 
length and almost negligible for 15 < x < 25 mm. Considering these results and the technological possibilities 
to build a SDD with a pitch smaller than 175 fim, we suggest 145 fim as the optimal pitch for the LOFT Wide 
Field Monitor. 

In order to characterize the detector spatial resolution with high precision we performed a new set of sim- 
ulations with a fixed anode pitch of 145 fim. In contrast to the simulations described above, this new set was 
performed by dividing each anode in 200 small sub-pixels, each with a surface of 14.5 /mix 1.75 mm, and generat- 
ing 10 3 photons of different energies (2.0, 2.5, 3.0, 4.5, 6.0, 10.0 and 20.0 keV) with a uniform spatial distribution 
inside the sub-pixel. In the left panel of Figure [7] we show the SDD spatial resolution FWHM (averaged on each 
sub-pixel) along the drift direction for 2 keV photons as a function of the absorption point in both the anodic and 
drift directions. In the right panel of the same figure we display the anodic spatial resolution. Such resolution is 
in general of the order of 20-30 fim and, as expected, becomes comparable with pitch/2 when the electron cloud 
is completely collected by one single anode. 

6. DISCUSSION AND CONCLUSIONS 



We developed a Monte Carlo simulator to study the effects of the charge diffusion on the spatial resolution of 
the SDDs. The simulated detector-FEE model takes into account all the most important factors which affect the 
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Figure 5. SDD spatial resolution (fim FWHM) for E p h=2 keV (upper panels) and E p h=3.5 keV (lower panels) for several 
photon absorption positions in the anodic and drift directions and for two different positions in the anodic direction (left 
panels: y=0, right panels: y=50 jj,m). The curves represent the resolution for 145 /J,m (black), 175 /im (green), 200 fim 
(brown), 225 /im (dark red), 250 /im (light red), and 294 fim (magenta) anode pitch. 
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Figure 6. SDD spatial resolution (/im FWHM) for E p h=6 keV for several photon absorption positions in the anodic and 
drift directions and for two different positions in the anodic direction (left panels: y=0, right panels: y=50 /im). The 
curves represent the resolution for 145 /im (black), 175 /im (green), 200 /im (brown), 225 /im (dark red), 250 (light 
red), and 294 (im (magenta) anode pitch. For energies E>6 keV, the anodic resolution becomes of the order of ~10 fim 
(x > 10 mm) for all the anode pitches considered in the simulations. 



detector imaging capabilities, investigating a number of configurations involving different detector designs and 
environmental conditions. The simulations have driven the design of the SDDs for the LOFT/WFM instrument, 
demonstrating that a pitch of 145 /jm is the optimal choice for maximizing the detector imaging performance. 
Such performance will be verified in the next months by using an ASIC-based FEE and the SDD prototypes 
already developed within the XDXL project in collaboration with Fondazione Bruno Kessler (FBK), Trento, 
Italy. With this design, extremely high spatial resolution (better than 70 /im FWHM in the anodic direction) 
can be obtained by fitting the charge cloud distribution with a discrete Gaussian function. The same procedure 
also allows to determine the photon absorption point in the drift direction with a resolution (FWHM) better 
than 8 mm at 2 keV, 5 mm at 3.5 keV and 3 mm at 6 keV. The Silicon Drift Detectors can thus be effectively 
used as a soft (E>2 keV) X-ray imaging detector, allowing to perform bi-dimensional photon-by-photon imaging 
in the energy range 2-50 keV, with low resource requirements in term of power, weight and complexity of the 
read-out electronics. 
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